Hadronic r decays and QCD 
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Abstract: Hadronic decays of the r lepton provide a clean source to study hadron dynamics in an 
energy regime dominated by resonances, with the interesting information captured in the spectral 
functions. Recent results on exclusive channels are reviewed. Inclusive spectral functions are the 
basis for QCD analyses, delivering an accurate determination of the strong coupling constant and 
quantitative information on nonperturbative contributions. Strange decays yield a determination of 
the strange quark mass. 



> 

o 

(N 

CD 



X 



1. Introduction 

Hadrons produced in r decays are borne out of 
the charged weak current, i.e. out of the QCD 
vacuum. This property garantees that hadronic 
physics factorizes in these processes which are 
then completely characterized for each decay chan- 
nel by spectral functions as far as the total rate is 
concerned . Furthermore, the produced hadronic 
systems have / = 1 and spin-parity = 0+, 1^ 
(V) and = Q- . 1+ (A). The spectral functions 
are directly related to the invariant mass spectra 
of the hadronic final states, normalized to their 
respective branching ratios and corrected for the 
T decay kinematics. For a given spin-1 vector 
decay, one has 



vi{s) = 
dNv 
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where Vud denotes the CKM weak mixing matrix 
element and S'ew = 1.0194 ± 0.0040 accounts for 
electroweak radiative corrections . 

Isospin symmetry (CVC) connects the r and 
e^e'" annihilation spectral functions, the latter 
being proportional to the R ratio. For example, j^^ve been pubhshedj|. New results from CLEO 



Radiative corrections violate CVC, as contained 
in the S'ew factor which is dominated by short- 
distance effects and thus expected to be essen- 
tially final-state independent. 

Hadronic r decays are then a clean probe of 
hadron dynamics in an interesting energy region 
dominated by resonances. However, perturbative 
QCD can be seriously considered due to the rel- 
atively large r mass. Many hadronic modes have 
been measured and studied, while some earlier 
discrepancies (before 1990) have been resolved 
with high-statistics and low-systematics exper- 
iments. Samples of ~ 4 x 10^ measured decays 
are available in each LEP experiment and CLEO. 
Conditions for low systematic uncertainties are 
particularly well met at LEP: measured samples 
have small non-r backgrounds (< 1%) and large 
selection efficiency (92%), for example in ALEPH. 

Recent results in the field are discussed in 
this report. 

2. Specific final states 
2.1 Vector states 

The decay r — > VtT^~'jt^ is now studied with large 
statistics of ~ 10^ events. Data from ALEPH 



are now available 



with the mass spectrum 



^i,x-{-^) (1-2) given in figure |] dominated by the p{770) res- 
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onance. Good agreement is observed between 
the ALEPH and CLEO data and the hneshape 
fits show strong evidence for the contribution of 
p(1400) through interference with the dominant 
amphtude. Fits also include a p(1700) contribu- 
tion, taken from e'^e~ data as the value of the t 
mass does not allow r data alone to tie down the 
corresponding resonance parameters. Thanks to 
the high precision of the data, fits are sensitive to 
the exact form of the p lineshape, with a prefer- 
ence given to the Gounaris-Sakurai parametriza- 
tion [Q over that of Kiihn-Santamaria ||] . 
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Figure 1: Mass distribution in CLEO r — > UTir~Tv° 
events. The solid curve overlaid is the result of the 
Kuhn-Santamaria fit, while the dashed curve has the 
p(1400) contribution turned off. 

The quality of data on e^e~ — > tt+tt^ has 
also recently improved with the release of the 
CMD-2 results from Novosibirsk |6|. A compar- 
ison of the mass spectrum as measured in e"'"e~ 
and r data is given in figure || (for this exercise 
the p — Lj interference has to be artificially intro- 
duced in the r data). Although the agreement 
looks impressive, it is possible to quantify it by 
computing a single number, integrating over the 
complete spectrum. It is convenient for this to 
use the branching ratio B{t Vr'iT~TT^) as di- 
rectly measured in t decays and computed from 
the e+e^ spectral function under the assump- 
tion of CVC. Using B{t ^ Vrh-TT°) ^ (25.79 ± 
0.15)% and subtracting out B{t ^ VrK^ii^) = 
(0.45 ± 0.02)% §, one gets B{t ly^T^-n^) = 
(25.34 ± 0.15)%, somewhat larger than the CVC 
value using all available e+e" data, -82^^'" — 



(24.65±0.35)% §. This 1.8cr discrepancy should 
be further investigated with a detailed examina- 
tion of the respective possible systematic effects, 
such as radiative corrections in e'^e~ data and 
7r° reconstruction in t data. CVC violations are 
of course expected at some point: hadronic vio- 
lation should be very small (~ {rriu — rudf' /m^), 
while significant effects could arise from long- 
distance radiative processes. Estimates show that 
the difference between the charged and neutral 
p widths should only be at the level of (0.3 ± 
0.4)% |0|. 
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Figure 2: Cross section for e^e~ ^ ■k'^-k~ com- 
pared to the ALEPH r data using CVC with p — lo 
interference built in (shaded band). 

The 47r final states have also been studied 
|l|. Tests of CVC are severely hampered by 
large deviations between different e+e" experi- 
ments which disagree well beyond their quoted 
systematic uncertainties. A new CLEO analy- 
sis studies the resonant structure in the 37r7r° 
channel which is shown to be dominated by wtt 
and aiTT contributions. The cjtt spectral func- 
tion shown in figure ^ is in good agreement with 
CMD-2 results and it is interpreted by a sum of 
yO-like amplitudes. The mass of the second state 
is however found at (1523 ± 10)MeV, in contrast 
with the value (1406±14)MeV from the fit of the 
27r spectral function. This point has to be clar- 
ified. Following a limit of 8.6% obtained earlier 
by ALEPH |l|, CLEO sets a new 95% CL hmit 
of 6.4% for the relative contribution of second- 
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class currents in the decay r — > v^-t^ lj from the 
hadronic angular decay distribution. 
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Figure 3: Fits of the CLEO spectral function for 
2.2 Axial- vector states 

The decay r Vr'in is the cleanest place to 
study axial- vector resonance structure. The spec- 
trum is dominated by the 1"*" ai state, known to 
decay essentially through pir. A comprehensive 
analysis of the ■k~2-k^ channel has been presented 
by CLEO. First, a model- independent determi- 
nation of the hadronic structure functions gave 
no evidence for non-axial- vector contributions (< 
17% at 90% CL) Second, a partial- wave am- 
plitude analysis was performed [p^ : while the 
dominant pi: mode was of course confirmed, it 
came as a surprize that an important contribu- 
tion (- 20%) from scalars (ct, /o(1470), /2(1270)) 
was found in the 2ti system. 

The ai ^ 7r^27r'^ lineshape is displayed in 
figure H where the opening of the K*K decay 
mode in the total ai width is clearly seen. The 
derived branching ratio, B{ai — > K*K) = (3.3 ± 
0.5)% is in good agreement with ALEPH results 
on the KKt: modes which were indeed shown 
(with the help of e+e~ data and CVC) to be 
axial-vector (ai) dominated with B{ai K*K) = 
(2.6 ± 0.3)% 1^. No conclusive evidence for a 
higher mass state {a'^) is found in this analysis. 

3. Inclusive spectral functions 

The T nonstrange spectral functions have been 
measured by ALEPH |, |l| and OPAL jl§ . The 
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Figure 4: (a) Svr mass spectrum from CLEO in the 
-K~2-K° channel with the lineshape from the resonance 
fit (zoom on the K* K threshold). In (b), the y^s- 
dependent a\ width is plotted with the different con- 
tributions considered. 

procedure requires a careful separation of vector 
(V) and axial- vector (A) states involving the re- 
construction of niulti-TT*' decays and the proper 
treatment of final states with a KK pair. The 
V and A spectral functions are given in figures ^ 
and 1^, respectively. They show a strong resonant 
behaviour, dominated by the lowest p and ai 
states, with a tendancy to converge at large mass 
toward a value near the parton model expecta- 
tion. Yet, the vector part stays clearly above 
while the axial-vector one lies below. Thus, the 
two spectral functions are clearly not 'asymp- 
totic' at the T mass scale. 

The V + A spectral function, shown in fig- 
ure 1^ has a clear converging pattern toward a 
value above the parton level as expected in QCD. 
In fact, it displays a textbook example of global 
duality, since the resonance-dominated low-mass 
region shows an oscillatory behaviour around the 
asymptotic QCD expectation, assumed to be valid 
in a local sense only for large masses. This fea- 
ture will be quantitatively discussed in the next 
section. 



3 



Michel Davier 




2 2.5 3 
s (GeV^) 

Figure 5: Inclusive nonstrange vector spectral func- 
tion from ALEPH and OPAL. The dashed line is the 
expectation from the naive parton model 
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Figure 6: Inclusive nonstrange axial-vector spectral 
function from ALEPH and OPAL. The dashed line 
is the expectation from the naive parton model. 



4. QCD analysis of nonstrange r de- 
cays 

4.1 Motivation 

The total hadronic t width, properly normalized 
to the known leptonic width, 

r(T- hadrons" Vr) , . 
-Kr = — ^rr^ — ^ — (4.1) 

should be well predicted by QCD as it is an inclu- 
sive observable. Compared to the similar quan- 
tity defined in e^e~ annihilation, it is even twice 
inclusive: not only all produced hadronic states 
at a given mass are summed over, but an inte- 
gration is performed over all the possible masses 
from m,r to m^. 
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Figure 7: Inclusive V + A nonstrange spectral 
function from ALEPH. The dashed line is the ex- 
pectation from the naive parton model, while the 
solid one is from massless perturbative QCD using 



0.120. 



This favourable situation could be spoiled by 
the fact that the scale is rather small, so that 
questions about the validity of a perturbative ap- 
proach can be raised. At least two levels are to be 
considered: the convergence of the perturbative 
expansion and the control of the nonperturbative 
contributions. Happy circumstances make these 
contributions indeed very small |l^, [l^ . 

4.2 Theoretical prediction for Rt 

The imaginary parts of the vector and axial- vector 
two-point correlation functions n|^'^'^^^^(s), with 
the spin J of the hadronic system, are propor- 
tional to the r hadronic spectral functions with 
corresponding quantum numbers. The non-strange 
ratio Rt can be written as an integral of these 
spectral functions over the invariant mass-squared 
s of the final state hadrons ||l^ : 

Rriso) = 127r5Ew / - ( 1 - - ) (4.2) 
J So V So/ 



1 + 2—1 Imn(i) (s + ie) + ImH^") (s + ie) 
So 



By Cauchy's theorem the imaginary part of XI^"^^ 
is proportional to the discontinuity across the 
positive real axis. 

The energy scale sq for sq — is large 
enough that contributions from nonperturbative 
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effects be small. It is therefore assumed that one 
can use the Operator Product Expansion (OPE) 
to organize perturbative and nonperturbative con- 
tributions ||2^ to Rt{sq). 

The theoretical prediction of the vector and 
axial- vector ratio Rry/A can thus be written as: 



Rry/A — -j\yud 



^S'ew 



(4.3) 



(D) 
d.yjA 



r>=4,6,8 



with the residual non-logarithmic electroweak cor- 



rection 5p 



EW 



lowing, and the dimension D 

c(2— mas 



0.0010 neglected in the fol- 
2 contribution 



(^i J ^fT^ from quark masses which is lower than 



0.1% for Ujd quarks. The term J'-*'-' is the purely 
perturbative contribution, while the 5^^^ are the 
OPE terms in powers of Sq of the following 
form 



^udyjA 



dimO=D 



(4.4) 



where the long-distance nonperturbative effects 
are absorbed into the vacuum expectation ele- 
ments iPud)- 

The perturbative expansion (FOPT) is known 
to third order j2^. A resummation of all known 
higher order logarithmic integrals improves the 
convergence of the perturbative series (contour- 
improved method FOPTci} ||2^. As some ambi- 
guity persists, the results are given as an average 
of the two methods with the difference taken as 
systematic uncertainty. 

4.3 Measurements 

The ratio is obtained from measurements of 
the leptonic branching ratios: 



Rr = 3.647 ±0.014 



(4.5) 



using a value B{t~ e~ Ve^r) = (17.794 ± 
0.045)% which includes the improvement in ac- 
curacy provided by the universality assumption 
of leptonic currents together with the measure- 
ments of B(t~ —f e~ C'e^^r), B(t~ — > /i~ V^Vr) 

and the r lifetime. The nonstrange part of Rr is 
obtained subtracting out the measured strange 
contribution (see last section). 



ALEPH 


Qs(m?) 


5np 


V 
A 


0.330 ±0.014 ±0.018 
0.339 ±0.013 ±0.018 


0.020 ± 0.004 
-0.027 ± 0.0004 


V+A 


0.334 ± 0.007 ±0.021 


-0.003 ± 0.0004 



Table 1: Fit results of Qs(mJ) and the OPE nonper- 
turbative contributions from vector, axial- vector and 
(V + A) combined fits using the corresponding ratios 
Rt and the spectral moments as input parameters. 
The second error is given for theoretical uncertainty. 

Two complete analyses of the V and A parts 
have been performed by ALEPH ||l5| and OPAL |l6 
Both use the world-average leptonic branching 
ratios, but their own measured spectral functions. 
The results on as{m'^) are therefore strongly cor- 
related and indeed agree when the same theoreti- 
cal prescriptions are used. Here only the ALEPH 
results are given. 

4.4 Results of the fits 

The results of the fits are given in table The 
limited number of observables and the strong cor- 
relations between the spectral moments intro- 
duce large correlations, especially between the 
fitted nonperturbative operators. 

One notices a remarkable agreement within 
statistical errors between the as(TO^) values us- 
ing vector and axial-vector data. The total non- 
perturbative power contribution to Rr.v+A is com- 
patible with zero within an uncertainty of 0.4%, 
that is much smaller than the error arising from 
the perturbative term. This cancellation of the 
nonperturbative terms increases the confidence 
on the as{mr) determination from the inclusive 
(V + A) observables. 

The final result is : 



a,{ml) = 0.334 ± 0.007cxp ± 0.021 



Lth 



(4.6) 



where the first error accounts for the experimen- 
tal uncertainty and the second gives the uncer- 
tainty of the theoretical prediction of Rt and the 
spectral moments as well as the ambiguity of the 
theoretical approaches employed. 

4.5 Test of the running of a.s{s) at low en- 
ergies 

Using the spectral functions, one can simulate 
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the physics of a hypothetical r lepton with a mass 



y/so smaller than m-r through equation (4.2) and 
hence further investigate QCD phenomena at low 
energies. Assuming quark-hadron duality, the 
evolution of Rt{so) provides a direct test of the 
running of Us (sq) , governed by the RGE /3- function. 
On the other hand, it is a test of the validity of 
the OPE approach in t decays. 
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Figure 8: The ratio Rt,v+a versus the square "r 
mass" So. The curves are plotted as error bands to 
emphasize their strong point-to-point correlations in 
So • Also shown is the theoretical prediction using the 
results of the fit at so — rri^. 
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Figure 9: The running of as (so) obtained from the 
fit of the theoretical prediction to _Rt,v+a(so). The 
shaded band shows the data including experimental 
errors. The curves give the four-loop RGE evolution 
for two and three flavours. 

The functional dependence of Rry+A{so) is 
plotted in figure || together with the theoretical 
prediction using the results of table |l]. Below 
1 GeV^ the error of the theoretical prediction 
of Rry+Aiso) starts to blow up due to the in- 
creasing uncertainty from the unknown fourth- 
order perturbative term. Figure |^ shows the plot 
corresponding to Fig. ||, translated into the run- 
ning of as(so), i.e., the experimental value for 
as(so) has been individually determined at ev- 
ery So from the comparison of data and theory. 



Good agreement is observed with the four-loop 
RGE evolution using three quark flavours. 

The experimental fact that the nonpertur- 
bative contributions cancel over the whole range 
1.2 GeV^ < So < TO^ leads to confidence that 
the as determination from the inclusive (V + A) 
data is robust. 

4.6 Discussion on the determination of 0^(771^) 

The evolution of the asim^) measurement from 
the inclusive {V + A) observables based on the 
Runge-Kutta integration of the differential equa- 
tion of the renormalization group to N^LO |23, 
|6| yields 

a,(Af|) = 0.1202±0.0008oxp±0.0024th±0.0010cvoi 

(4.7) 

where the last error stands for possible ambigu- 
ities in the evolution due to uncertainties in the 
matching scales of the quark thresholds Pq| . 




Energy (GeV) 

Figure 10: Evolution of the strong coupling mea- 
sured at m^) to Mz predicted by QCD compared 
to the direct measurement. The evolution is carried 
out at 4 loops, while the flavour matching is accom- 
plished at 3 loops at 2 rric and 2 nit thresholds. 



The result (4.7) can be compared to the de- 
termination from the global electroweak fit. The 
variable Rz has similar advantages as Rr, but 
it differs concerning the convergence of the per- 
turbative expansion because of the much larger 
scale. It turns out that this determination is 
dominated by experimental errors with very small 
theoretical uncertainties, i.e. the reverse of the 
situation encountered in r decays. The most re- 
cent value|3 yields as(M|) = 0.119 ± 0.003, 
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in excellent agreement with ([l.7|). Fig. |l^ illus- 
trates well the agreement between the evolution 
of as{m'^) predicted by QCD and as(Af|). 

5. Applications to hadronic vacuum 
polarization 

5.1 Improvements to the standard calcula- 
tions 

Prom the studies presented above we have learned 
that: 

• the 7=1 vector spectral function from t 
decays agrees with that from e~^e~ annihi- 
lation, while it is more precise for masses 
less than 1.6 GeV as can be seen on fig- 
ure Small CVC violations are expected 
at a few 10"'^ level from radiative p 
decays and SU(2)-breaking in the n and p 
masses. 

• the description of Rr by perturbative QCD 
works down to a scale of 1 GeV. Nonper- 
turbative contributions at 1.8 GeV are well 
below 1 % in this case. They are larger (~ 3 
%) for the vector part alone, but reason- 
ably well described by OPE. The complete 
(perturbative -I- nonperturbative) descrip- 
tion is accurate at the 1 % level at 1.8 GeV 
for integrals over the vector spectral func- 
tion such as i?r,y. 



ii| . x"^(V, I=l)v, (ALEPH) - 

h + - - 

1 ~ e e -»(V, 1=1) 
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Figure 11: Global test of CVC using r and e'^e 
vector spectral functions. 



These two facts have direct applications to 
calculations of hadronic vacuum polarisation which 
involve the knowledge of the vector spectral func- 
tion: the muon magnetic anomaly and the run- 
ning of a. In both cases, the standard method in- 
volves a dispersion integral over the vector spec- 
tral function taken from the e~^e~ — > hadrons 
data. Eventually at large energies, QCD is used 
to replace experimental data. Hence the preci- 
sion of the calculation is given by the accuracy 
of the data, which is poor above 1.5 GeV. Even 
at low energies, the precision can be significantly 
improved at low masses by using r data [lO[ |. 

The next breakthrough comes about using 
the prediction of perturbative QCD far above 
quark thresholds, but at low enough energies (com- 
patible with the remarks above) in place of non- 
competitive experimental datapsf. This proce- 
dure involves a proper treatment of the quark 
masses m the QCD prediction[||. 

Finally, it is still possible to improve the con- 
tributions from data by using analyticity and 
QCD sum rules, basically without any additional 
assumption. This idea, advocated in Ref. p9| , 
has been used within the procedure described 
above to still improve the calculations Q. 

The experimental results of R{s) and the the- 
oretical prediction are shown in figure The 
shaded bands depict the regions where data are 
used instead of theory to evaluate the respec- 
tive integrals. Good agreement between data and 
QCD is found above 8 GeV, while at lower ener- 
gies systematic deviations are observed. The R 
measurements in this region are essentially pro- 
vided by the 772 ^ and MARK I || collabo- 
rations. MARK I data above 5 GeV lie system- 
atically above the measurements of the Crystal 
Ball III and MDl || Collaborations as weU as 
the QCD prediction. 

5.2 Muon magnetic anomaly 

By virtue of the analyticity of the vacuum polar- 
ization correlator, the contribution of the hadronic 
vacuum polarization to can be calculated via 
the dispersion integral p5[ 

00 

= i / ds<J^adis)K{s) (5.1) 

4mJ 
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Aa(s)==-47rQ!Re [U'^{s) - U'^{0)], one has 
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Figure 12: Inclusive hadronic cross section ratio 
in e"'"e~ annihilation versus the cm. energy -^5. 
Additionally shown is the QCD prediction of the 
continuum contribution from Ref. j28| as explained 
in the text. The shaded areas depict regions were 
experimental data are used for the evaluation of 
Aahad(M|) and ajj'"* in addition to the measured 
narrow resonance parameters. The exclusive e"'"e~ 
cross section measurements at low cm. energies are 
taken from DMl,DM2,M2N,M3N,OLYA,CMD,ND 
and r data from ALEPH (see Ref. for detailed 
information) . 



Here (Thad(s) is the total e+e~— > hadrons cross 
section as a function of the cm. energy-squared 
s, and K{s) denotes a well-known QED kernel. 

The function K{s) decreases monotonically 
with increasing s. It gives a strong weight to the 



low energy part of the integral (5.1). About 91% 
of the total contribution to al^^*^ is accumulated 
at cm. energies y^s below 2.1 GeV while 72% of 
a^^'^ is covered by the two-pion final state which 
is dominated by the p(770) resonance. The new 
information provided by the ALEPH 2- and 4- 
pion spectral functions can significantly improve 
the a^.'^'^ determination. 

5.3 Running of the electromagnetic cou- 
pling 

In the same spirit we evaluate the hadronic con- 
tribution Aa(s) to the renormalized vacuum po- 
larization function li'yis) which governs the run- 
ning of the electromagnetic coupling a{s). With 



a{0) 



1 - Aa{s) 



(5.2) 



where 47ra(0) is the square of the electron charge 
in the long- wavelength Thomson limit. 

The leading order leptonic contribution is 
equal to 314.2 x 10"**. Using analyticity and uni- 
tarity, the dispersion integral for the contribution 
from the light quark hadronic vacuum polariza- 
tion Aahad(M|) reads pi 



47r2 Q 



Re 



ds ■ 



4mi 



-Mi- 



(5.3) 

where a{s) — 167r^a^ (s)/s • ImH^(s) from the 
optical theorem. In contrast to a^'^'^, the inte- 
gration kernel favours cross sections at higher 
masses. Hence, the improvement when includ- 
ing T data is expected to be small. 

The top quark contribution can be calculated 
using the next-to-next-to-leading order pre- 
diction of the total inclusive cross section ratio 
R from perturbative QCD |2^, The evalua- 
tion of the integral ( |5.3| ) with mtop — 175 GeV 
yields Aatop(A^I) = -0.6 x IQ-*. 

5.4 Results 

The combination of the theoretical and experi- 



mental evaluations of the integrals (5_^) and (5.1) 
yields the results 

Aahad(Afl) - (276.3 ± l.lexpi 1.1th) x 10"^ 
a-i(A'/|) = 128.933 ±0.015cxp± 0.015th 
a"^;^ = (692.4 ± 5.6exp ± 2.6th) x 10"!° 
a™ = (11659159.6 

±5.6exp ± 3.7th) X 10-1° (5.4) 



and a^""^ = (187.5±1.7oxp±0.7th) x IQ-i"^ for the 
leading order hadronic contribution to a^. The 
total value includes additional contributions 
from non-leading order hadronic vacuum polar- 
ization (summarized in Refs.[^ p^ ) and light- 
by-hght scattering ||3^, Q contributions. Fig- 
ures |l^ and ^ show a compilation of published 
results for the hadronic contributions to a(M|) 
and a^. Some authors give the hadronic contri- 
bution for the five light quarks only and add the 
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top quark part separately, 
rected for in Fig. fsl 



This has been cor- with the 'standard' a(M|) and 

log Mh - 1.97l°:i (5-6) 
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Figure 13: Comparison of Aahad(A/|) evaluations. 
The values are taken from Refs. ^ 0, 

[10[§, §, §,§. 
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Figure 14: Comparison of a^"^ evaluations. The 
values are taken from Refs. ||, ^ 1^, 0, 



5.5 Outlook 

These results have direct implications for phe- 
nomenology and on-going experimental programs. 

Most of the sensitivity to the Higgs boson 
mass originates from the measurements of asym- 
metries in the process e'^e~ — > Z — > fermion 
pairs and in fine from {si'n?Ow)ef f — ■ Un- 
fortunately, this approach is limited by the fact 
that the intrinsic uncertainty on a{M^) in the 
standard evaluation is at the same level as the 
experimental accuracy on s^. The situation has 
completely changed with the new determination 
of a(M|) which does not limit anymore the ad- 
justment of the Higgs mass from accurate ex- 
perimental determinations of sin^^w The im- 
provement in precision can be directly appreci- 
ated on the relevant variable log Mh with Mh 
in GeV/c2 p^: 



log Mh 



1. 



j+0.31 
'-0.39 



(5.5) 



with the QCD-improved value |3(][ |. 

The interest in reducing the uncertainty in 
the hadronic contribution to a^'^'^ is directly linked 
to the possibility of measuring the weak contri- 
bution: 



(5.7) 



10 



where a^^^ ^ (11658 470.6 ± 0.2) x 10" 
the pure electromagnetic contribution (see pl| 
and references therein), a^'^ is the contribution 
from hadronic vacuum polarization, and a^"^""^ = 
(15.1 ± 0.4) X 10-^° (H, H, H accounts for 
corrections due to the exchange of the weak in- 
teracting bosons up to two loops. The present 
value from the combined /i^ and measure- 



ments 54 



= (11 659 230 ±85) x 10 



-10 



(5.8) 



should be improved to a precision of at least 
4 X 10^^" by a forthcoming Brookhaven exper- 
iment (BNL-E821) well below the expected 
weak contribution. Such a programme makes 
sense only if the uncertainty on the hadronic term 
is made sufficiently small. The improvements de- 
scribed above represent a significant step in this 
direction. 

6. Strange r decays and rUg 

6.1 The strange hadronic decay ratio i?T,s 

As previously demonstrated in Ref. [ p6| , the in- 
clusive r decay ratio into strange hadronic final 
states, 



hadrons 



s=-i 



r(r- 



(6.1) 



can be used due to its precise theoretical pre- 
diction ||l^, ^ to determine ms(s) at the scale 
s = M^. Since then it was shown that 
the perturbative expansion used for the massive 
term in Ref. ||5^ was incorrect. After correction 
the series shows a problematic convergence be- 
haviour p8, p£ 
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Figure 15: Some of the new ALEPH results on 
decay modes with kaons compared to the published 
data. 



Similarly to the nonstrange case, the QCD 



prediction is given by equations (4^) and (4.3) 
where the attention is now turned to the ,5(2-mass) 
term, important for the relatively heavy strange 
quark. The corresponding perturbative expan- 
sion is known to second order for the J = 1 + 
part and to the third order for J = |5^, |5^ . 
While the J = 1+0 series behaves well, the J = 
expansion in fact diverges after the second term. 

Following these observations, two methods 
can be considered in order to determine ms(m^): 



• in the inclusive method, the inclusive strange 
hadronic rate is considered and both J = 
1 + and J = are included with their re- 
spective convergence behaviour taken into 
account in the theoretical uncertainties. 



6.1.1 Nevif ALEPH results on strange de- 
cays 

ALEPH has recently published a comprehensive 
study of r decay modes including kaons (charged, 

and if^) H H, I up to four hadrons in 
the final state. A comparison with the published 
results is given in figure ^ 

The total branching ratio for t into strange 
final states, Bs, is 

= (2.87+0.12)% (6.2) 

corresponding to 

i?^ 5 = 0.1610+0.0066 (6.3) 
Since the QCD expectation for a massless quarks 



is 0.1809 + 0.0036, the result (|6j) is evidence for 
a massive s quark. 

The strange spectral function is given in fig- 
ure |l^, dominated at low mass by the K* reso- 
nance. 



ALEPH 

□ K7t 

m kkk 

m K37i+K"Ti (MC) 
m K4n (MC) 

□ K57t (MC) 




Mass^ (GeW/cY 



Figure 16: The strange spectral function measured 
by ALEPH. 



the '1+0' method singles out the well-behaved 
J = 1 + part by subtracting the exper- 
imentally determined J = longitudinal 
component from data. The measurement 
is then less inclusive and the sensitivity to 
uis is significantly reduced; however, the 
62 perturbative expansion is under control 
and the corresponding theoretical uncer- 
tainty is reduced. 



6.1.2 The ALEPH analysis for rus 

As proposed in Ref. Q and successfully applied 
in cks(m^) analyses, the spectral function is used 
to construct moments 



T)kl — 



(6.4) 
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In order to reduce the theoretical uncertain- 
ties one considers the difference between non- 
strange and strange spectral moments, properly 
normalized with their respective CKM matrix el- 
ements: 



1 



nkl 

IT/ \2^r,S=a 



1 



\Vus 



T)kl 



(6.5) 



for which the massless perturbative contribution 
vanishes so that the theoretical prediction now 

reads (setting m„ = rud = 0) 



A^' = as*, 



EW 



-S' 



fc^ (2 — mass) 



- E 

0=4,6, 



~gkl(D) 



(6.6) 



For the CKM matrix elements the values \ Vud\ = 
0.9751 ± 0.0004 and \Vu,\ = 0.2218 ± 0.0016 H 
are used, while the errors are included in the the- 
oretical uncertainties. 

Figure |l^ shows the weighted integrand of 
the lowest moment A^*' from the ALEPH data, as 
a function of the invariant mass-squared, and for 
which the expectation from perturbative QCD 
vanishes. 
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Figure 17: Integrand of Eq. (|6j) for (k=0, 1=0), 
i.e., difference of the Cabibbo corrected non-strange 
and strange invariant mass spectra. 

Subtracting out the experimental J ~ con- 
tribution (essentially given by the single K chan- 
nel), a fit to 5 moments is performed with the 
safer J = 1 -I- method, yielding 



ms(m^) 



(176 



h37o 

-48o: 



:^28:::) Mev/c^ 



(6.7) 



6) 



0.039 ± 0.016exp ± 0.014th (6.8) 



^(8) ^ -0.021 ± O.OMexp ± 0.008th(6.9) 

with additional (smaller) uncertainties from the 
fitting procedure and the determination of the 
J — part. The quoted theoretical errors are 
mostly from the Vus uncertainty. The D ~ Q 
and D = 8 strange contributions are found to be 
surprizingly larger than their nonstrange coun- 
terparts. If the fully inclusive method is used 
instead (with the problematic convergence) the 
result msiml) = {UdtltZl-llTJ ^eV/c^ is ob- 
tained. 



The result ( |6.7D can be evolved to the scale of 
1 GeV using the four- loop RGE 7-function |66| , 
yielding 



m,(l GeV^) = (234 i^^) MeV/c^ 



(6.10) 



This value of rn^ is somewhat larger than 
previous determinations pTf , but consistent with 
them within errors. 

7. Conclusions 

The decays t ^ + hadrons constitute a clean 
and powerful way to study hadronic physics up 
to i/i ~ 1.8 GeV. Beautiful resonance analyses 
have already been done, providing new insight 
into hadron dynamics. Probably the major sur- 
prize has been the fact that inclusive hadron pro- 
duction is well described by perturbative QCD 
with very small nonperturbative components at 
the T mass. Despite the fact that this low-energy 
region is dominated by resonance physics, meth- 
ods based on global quark-hadron duality work 
indeed very well. 

The measurement of the vector and axial- 
vector spectral functions has provided the way 
for quantitative analyses. Precise determinations 
of as agree for both spectral functions and they 
also agree with all the other determinations from 
the Z width, the rate of Z to jets and deep in- 
elastic lepton scattering. Indeed from r decays 



as{Ml)r = 0.1202 ±0.0027 



(7.1) 



in excellent agreement with the average from all 
other determinations [BSl 



Q^s(-^^z)nori-r 



0.1187 ±0.0020 



(7.2) 
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The use of the vector r spectral function 
and the QCD-based approach as tested in r de- 
cays improve the calculations of hadronic vac- 
uum polarization considerably. Significant re- 
sults have been obtained for the running of a 
to the Z mass and the muon anomalous mag- 
netic moment. Both of these quantities must be 
known with high precision as they give access to 
new physics. 

Finally the strange spectral function has been 
measured, providing a determination of the strange 
quark mass. 
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